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This equation of state may be useful for interpolation and for computations which include thermo- 
dynamic properties of parahydrogen. A principal advantage is that it employs only one set of 24 
coefficients for all fluid states. The basic form is relatively simple. An adjustment term for the critical 
region leads to good agreement with published PVT data and specific heats in this especially difficult 
area. The forms of fitting functions were developed by systematic trial methods using least squares 
to find the coefficients. 
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List of Symbols and Units 

A(p), B(p), C{p) density-dependent parameters. 

CY, specific heat at constant volume; C° for ideal gas 

states. 
<i, density. 
d c , d t , critical-point, and triple-point densities. 

A„s(l00/n) •Sl^/Pca.c-ll for n data. 

i 

J, the joule. 

/, the liter, defined as 1000 cm ;i (12th Gen. Conf. Wts. 

Measures, Paris, October 1964).** 
P, pressure. 

P ( , Pt, critical-point, and triple-point pressures. 
Pt = 0.0695 atm [6|. 

Il(p), a density-dependent, parametric pressure. 
R, the gas constant, 82.0597 cm'* atm/g mol deg[3|. 
p = d/d^ density reduced at the triple point. 
a = djdc, density reduced at the critical point. 
7\ temperature, deg Kelvin (NBS 1955 low-temp, scale). 
TV, T ( critical-point, and triple-point temperatures. 
t = 27 7V, temperature reduced at the critical point. 
v = 1/rf, molal volume. 

#=s 77 7 1 /, temperature reduced at the triple point. 
y, a defined function of P, II, 7\ d. 
Z(cr, t), a defined function. 
(pier, r), a defined function. 



*Cryogenics Division. NHS Laboratories, National Bun 
80302. 
**NBS Misc. Publ. 279, p. 203 (1965). 
1 Figures in brackets indicate the literature referenc 



of Standards, Boulder. Colo 



i at the end of llii* paper. 



1. Introduction 

Our initial computations of thermodynamic proper- 
ties of parahydrogen fl| l were accomplished hy use of 
the Benedict-Wehh-Rubin equation of state as modified 
by Strobridge [2|. To obtain the required accuracy, 
however, it was necessary to evaluate two sets of con- 
stants for two regions <>i the P—V—T surface. Serious 
computational complications then arose from the need 
to establish a practical continuity of thermodynamic 
properties over the boundaries of these regions. 
Partly for this reason, we used polynomial methods in 
our formal publication on the thermodynamic proper- 
ties [3|. Despite these tabulations, there remains a 
need for an equation of state which can be used for 
interpolation, and for a multitude of computational 
investigations. A prerequisite of the following work 
has been that the equation of state shall require but a 
single set of constants for all fluid states of the sub- 
stance (to avoid the above-mentioned problem of 
continuities). 

Well-behaved descriptions of P—p—T behavior are 
obtained with the form P(p, T), that is, with density 
and temperature as arguments or parameters. Den- 
sity-dependence of specific heats then is obtained from 
P—p-T data via the thermodynamic relation. 

(dC v ldp) T =(-Tlp*)(d*Pldr%. 

Thus an exacting requirement upon any equation of 
state is to give the curvature of isochores, (d 2 PldT~) p , 
accurately as a function of density and temperature. 
We have elected to develop an isochoric form of 
equation, i.e., an expression for pressure in terms of 
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temperature. As any accurate equation for high 
densities must be empirical, we considered that tem- 
perature dependence along isochores (fixed mean 
intermolecular distance) probably would be less com- 
plicated than density-dependence along isotherms [4], 
We have allowed the density-dependence of parameters 
in the equation to be given by P—p—T data on isochores. 
These data are interpolated from experimental 
isotherms [5J. 

2. The Behavior of P — p — T Data on 
Isochores 

In this section we survey the set of interpolated data 
on isochores with a view to discovering any useful 
generalities. This survey is possible because the 
experimental data on isotherms are dense, covering a 
wide range of conditions, and because they have been 
interpolated quite accurately on isochores [5]. For 



convenience in obtaining derivatives for this survey, 
we take values computed from polynomial representa- 
tions of the isochores with coefficients given by [3J. 
The derivatives therefore may be of low accuracy in 
some regions. 

Plots of P versus T at constant density (isochores) 
appear to be a family of straight lines. (This empha- 
sizes the need for extreme accuracy in P — p — T data 
used to obtain derivatives.) First derivatives (dP/dTjp, 
however, reveal useful characteristics. Figure 1 
gives dimensionless slopes, v • (dP/dT)/R, versus T on 
isochores. For d > 0.013 g mol/cm 3 each plot has a 
maximum, corresponding to an inflection in the P(T) 
relation. Curvatures, (d 2 P/dT J2 )p, are positive at low 
temperatures and negative at high temperatures. 
The locus of these inflections is important for the 
computation of specific heats on isotherms, due to 
the inversion of sign of the integrand. Figures 2 and 
3 give the approximate P—T and T—p coordinates of 
the locus. On isotherms from 33 to about 52 °K we 
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FIGURE 1. Slopes of selected isochores versus T. 
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FIGURE 2. Locus of isochore inflection points in P — T coordinates. 
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FIGURE 3. Locus of isochore inflection points in T-p coordinates. 

The clashed line from left to right represents saturated vapor, saturated liquid, and freezing liquid, respectively. 
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see that {d 2 Plc)T 2 ) p might have three inversions of 
sign as the density increases from zero to very high 
values. Part of this behavior already is known [16|. 

A further important characteristic of isochores is 
seen only by careful examination of differences in the 
data. The terminating ranges near saturation (low 
temperatures) exhibit sharp curvatures {d 2 P/dT 2 ) p as 
T^> r sat , appropriately named "hooks." In the criti- 
cal region these curvatures are negative for p < p c 
and positive for p > p r . 

3. A Preliminary Equation 

We at first sought an elementary form for the sig- 
moid behavior of isochores. It should reduce to 
perfect gas behavior at low densities and also at high 
temperatures. We may include the purely density- 
dependent, parametric pressure found in the van der 
Waals and Beattie-Bridgeman equations, using the 
symbol U(p) [4J. 

Assume that the inflection-point temperature, 0(p), 
is known for each isochore. For each, define a new 
temperature scale, t(p,T) = T/0, and assume that 
0(0) = 0, (fig. 3). The simple form, 

[P+II(p)]/(J?-0-d)=*-exp [(l-i»)/(i» •*"■)], 

with 0< m # 1, then yields a derivative, (d 2 PldT 2 ) p ^ 
which is positive for T < and is negative for T > 0. 

Woolley et al., used an exponential form for hydro- 
gen at high temperatures, omitting the term Il(p) [6|. 
We have, therefore, investigated the general form 



v 

I 



log, Y=j?A r T~ ', 



(1) 



using the definition, 

Y = (P + Il)/(R-T'd). 

One must discover necessary values for exponents, 
Hi, and the corresponding number of terms, N. Den- 
sity-dependence of the A\ then is obtained from P(T) 
data by least squares under an iterative search for 
ll(p) on each isochore. 

Lengthy investigations failed to yield a behavior 
for ri(p) which was sufficiently regular to permit ana- 
lytical representation. With more than three terms 
(of the numerous powers in T investigated), the be- 
havior of ll(p) and all Ai{p) was exceedingly irregular. 
With three terms, the behavior of Il(p) above 0.01 
g mol/cm 3 was in rough accord with the Beattie- 
Bridgeman description, which we adopted for the con- 
strained form of II, 



n( P )-p 1 -p 2 -p 2 -p 3 . 



(2) 



When Pi and P 2 have been found (by iterative meth- 
ods), the A\ have a regular dependence on p. The 



lowest deviation, A p (see list of symbols), was obtained 
with the equation 



log, Y=A- p • x- ,/4 + B • p • x~ l + C • x~'\ 



(3) 



where x = TjT t , and symbols A, B, C, replace the A\ 
of eq (1). 

In each of the first two terms on the right one power 
of p has been introduced explicitly to diminish density- 
dependence of A and B. We note that exponent 
n\ = 1/4 is required for description of virial coefficients 
at low densities [7], and exponent 7*3 = 3 is found in 
the Beattie-Bridgeman equation. At sufficiently low 
densities, the term with ri^ — l combines with F[(p), 
independent of temperature. With Pi = 703.2813 atm 
and Pi = 88.3125 atm, one obtains the behavior for 
A(p), P(p), C{p) shown by figure 4. 
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Data needed for developing analytical descriptions 
of the parameters A(p), B(p), C{p) were obtained as de- 
scribed in appendix 6.1. The form of each descrip- 
tion, given below, is the result of a considerable in- 
vestigation, imposed by the need for high percision 
(appendix 6.2). In each, the number of terms and the 
values of exponents were explored systematically. In 
brief, the behaviors of figure 4 have been described 
as follows: 

A(p) =A X +A 2 'p *exp (— A 3 -p), (3-a) 

B(p)=B l + exp (B 2 + B 3 -p 

+ B 4 -p 2 + J B 5 -p 4 + ^-p 5 ), (3-b) 

C{p)=C x p + C 2 -p 2 + C z -p* + C 4 'p 4 

-hC 5 -p 5 + C (i -p 7 + C 7 -p 1() + C 8 -p 14 . (3-c) 

We found constants for eqs (3-a) to (3-c), omitting 
one datum at the lowest temperature (near saturation) 
from each of 85 isochores of interpolated data, 1.0 



^Sd^43.0 g mol/liter [5J. Mean deviation lor the 
remaining 2002 points was A,, = 0.026 percent. Ex- 
cessive deviations near the critical point, however, 
precluded calculation of specific heats in this region. 

Anticipating an adjustment function for the critical 
region, we redetermined a set of 19 constants for eq (3), 
excluding an additional 59 data in the domain, 32 
^ T ^ 35 °K, 8.5 ^ d ^ 22.5 g mol/liter. The pressure 
was constrained at the liquid triple-point by means of 
a relatively small adjustment to the value of con- 
stant Cx. 

These 19 constants, f\, . . ., Cs, are given in 
table 1. Mean deviations for eq (3) appear in table 2, 
and individual deviations (P — P C a\c) locale appear in 
table 3, given in parts per ten-thousand at tempera- 
tures through 42 °K. For higher temperatures they 
do not exceed a few parts in ten-thousand. (Where 
data do not exist, as in two-phase regions, or at P>350 
atm, tables 3 and 5 contain zero entries.) Figure 5 
gives the difference between left and right sides of 
eq (3), multiplied by 10 4 , on the 33-deg isotherm. 



Table 1. Constants for eqs (3) and (4) 



Equation (3) 


Equation (4) 


Constant 


Value 


Constant 


Value 


Constant 


Value 


T r °K 

T c °K 

d r g mol/cm 3 

P, atm 

Pi atm 

A x 


13.803 [6] 
32.976 [12] 
(1/26.176) [6] 

697.0 
81.0 
2.1350 2896 
1.0170 4119 
1.4268 4000 

-1.1427 1000 
1.9682 3016 

-2.2337 9180 


Bx 
B, 
Be 

Cx 

Ci 

c, 
c 4 

c, 

Ce 

c 7 

c H 


-0.2410 0077 

1.8319 9451 

-0.9445 7676 

-2.0176 0851 

-13.9844 7635 

128.3251 2808 

-294.2851 8710 

239.1130 5168 

-67.1458 8122 

10.7210 8689 

-0.7702 5591 


n 


8 


d r . ^ mol/cm 3 
C K 
Dx 
D-> 

Dx 
Ex 


0.0154 40 

-.7702 5578 

.0004 5110 

.0166 2606 

-.0105 4950 


At 


7.28 


A-, 


3.32 


Bx 




B» 




£, 









Table 2. Mean deviations for eqs (3) and (4) 



Equation number 


(3) 

19 

a 1943 


(4) 


Number of constants 


19 + 6 


No. of points used 


b 781 
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No. of points for deviation 


Mean deviation, A /; , % 
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a One point omitted at lowest temperature of each of 85 isochores, 
and all points in domain, 32 ^ T^ 35 °K, 8.5 ^ d ^ 22.5 g mol/liter. 

b For the six constants in Z{cr. r), data used in the domain. 17 ^ 7 1 
*S 40 °K, 1.0 ^ d *S 36.0 g mol/liter. 

( Points in the arbitrary critical domain, 32^7^35 °K, 8.5 =s d FIGURE 5. Difference between left and right sides of eq (3), multi- 
=S 22.5 g mol/liter. plied by 10 4 , on 33-deg isotherm. 
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4. The Adjusted Equation 

The deviation plot (fig. 5) shows behavior of an ad- 
justment on isotherms which must be made to the pre- 
liminary equation in the critical region. We see that 
it may correspond at least to an incipient level spot or 
shelf in the P(p, T) surface in three-dimensions, and 
we should not expect the preliminary equation to de- 
scribe this sharp irregularity. We have, therefore, 
added an adjustment term, Z(cr, r), a function of re- 
duced variables cr = d/d c , r = T/T c , as follows, 



loge Y=A • p-x^ + B' p'X- x + C'X-* + Z((j, r). 



(4) 



The adjustment diminishes rapidly at densities away 
from the critical point, and at temperatures above the 
critical point. 

Recent work suggests nonanalytic PVT behavior 
at the critical point, corresponding to a logarithmic 
infinity in C v [8, 9, 10, 11 J. The form of adjustment 
needed here, however, already is restricted by the pre- 
liminary equation of state. We have chosen an ana- 
lytic form, with the justification that it describes the 
experimental facts quite well in the range investigated. 

The form of Z, below, permits representation of 
the deviations in table 3 and also yields an acceptable 
set of specific heats in the critical region. With h\e 
coefficients D\, . . . , E\ the function is 



Z(cr, r) = D(or) • exp [— $(cr, t)] , 

D(a) =Z), -o-+(Z) 2 -o--hZ), -o- 2 )-(o--l), 



<A(o\t) = Z) 4 - (o— 1) 2 + £i -(t"-1) 



(4-a) 



(4-b) 
(4-c) 



We discuss briefly the basis for eq (4-a). Devia- 
tions at 33 deg (very near the critical isotherm at 
32.976 deg) require at least a cubic polynomial in 
density, with roots at o _ = 0, <x~l, cr > I: yet these 
deviations diminish rapidly at the highest densities. 
Equation (4-a) describes this behavior on isotherms. 
Temperature-dependence of these deviations is strong: 
when t=1.1 the deviations have about vanished. For 
this reason we introduced the factor exp(— E\ • r n ). It 
gives the sharp isochore curvatures or "hooks" which 
contribute to calculated specific heats in the critical 
region. It also allows the calculated deviations to 
persist at temperatures below critical at densities 
near the saturation envelope. 

On the right of eq (4-c) we selected the form of the 
first term, Z) 4 -(cr — l) 2 , following some tedious work 
with the form D 4 • (cr™ — 1 ) . Using the latter atm ~ 4, 
it was necessary to increase the powers of cr in eq 
(4-b) by two or more, e.g., 



D( 



Z), • O 3 + (D 2 ' cr 4 + D : , ■ cr 5 ) • (or - 1), 



signs of coefficients D\, D 2 , D$ exhibited radical be- 
haviors, and we abandoned this search. 

The constants of eq (4-a) have been determined for 
781 data in the domain 17 ^ T ^ 40 °K; 1.0 ^ d ^ 36.0 
g mol/liter, as described in section 5, below. Our 
general method was the following: We assumed values 
for exponent n, and then found Z) 4 , E\ by iterative 
procedures which minimize A ;> . In every step we 
used D\ for constraint to the assigned critical-point 
pressure, and /)•_>, D% to ensure that dP/dcr = d 2 P/dcr 2 =0 
at the assigned critical point. (An estimate of n ~ 10 
at first was obtained by assuming that exp [— Ei'(t h — 1)] 
~ 0.02 at t=1.1, and that its maximum positive curva- 
ture shall occur near t= 1.) 

We have used the critical temperature established 
by Roder et al., T c = 32.976 °K [12]. The critical pres- 
sure is computed at this temperature by means of 
the vapor-pressure equation of Weber et al., for ac- 
curate consistency, P ( = 12.759 263 atm [13j. The 
critical density has been varied within limits of un- 
certainty [14] to obtain consistency between specific 
heats from eq (4) and those tabulated in [3], as de- 
scribed below. 

5. Examination of the Equation 

An exacting test of the equation is behavior of calcu- 
lated specific heats, C r = C° v +C*, where the "residual 
specific heat" is 



C* - const. - (- D • (" (d 2 Pldr) p dp/p\ 
Jo 



(5) 



to obtain a deviation, A ;> , comparable with that found 
via eq (4-b) above. Even so the magnitudes and the 



with const. = (0.101325 J/cm :i a.lm)/{d ( ) when d t is in 
g mol/cm 3 . As a rough guide, deviations within 5 to 
10 percent of C v often are considered satisfactory. 
In this work we use the tabulations of [3] for compari- 
son, as they are in good agreement with experimental 
data. 

We have varied parameters n and d c of eq (4), 
examining behavior of the mean deviation A p and also 
of the mean deviation AC V between 470 specific heats 
from eq (4) and from [3]. Table 4(A) presents these 
results. In general, a minimum in A p does not coin- 
cide with the best minimum in AC r . We have made 
an arbitrary selection of n = 8, d c = 15.44 g mol/liter. 
Corresponding values for C 8 , D\, D>, Dz, Z) 4 , and E\ are 
given in table 1. 

Table 5 presents residual specific heats from eq (4) 
with above constants and compares them with values 
from [3]. The first row at d=0 gives C% from [6]. At 
T<T C the reference states for liquid are saturated 
liquid (not ideal gas), with properties given by table 6. 
Mean deviation for the 470 points of table 5 is 0.077 
J/g mol deg, or roughly half of 1 percent. The greatest 
deviations are about 3 percent in neighborhood of the 
critical point. 

It is interesting to see that integration on the 33 °K 
isotherm (within 0.02 °C of the critical point) yields 
satisfactory specific heats at the higher densities. 
Nevertheless, the maximum in C v on this isotherm 
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Table 3. Individual relative deviations from eq (3), 



DEN 
1.0 
1.5 
2.U 
2.5 
3.0 
3.b 
4.0 
4.b 
5.0 
5.5 
6.0 
6.b 
7.0 
7,b 
B.O 
B.b 

9.0 

9.5 
1".0 
10.5 
11.0 
11.5 
12.0 
12.5 
1J.0 
13.5 
14.0 
14.5 
15. 
15. b 
lfc. 
16.5 
17.0 
17. 5 
1H.0 
IB. 5 
19.0 
1<*.5 
20. 
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Table 4. Effects of parameters in eq [4) 
Part A, T c = 32.976 °K 



Table 4. 



Effects of parameters in eq (4) (continued) 
PartB. rc = 8 





d c , 
g mol/1 


(A p , %) • 100 


AC,, 
J/mol deg 




88 a 


781 b 


2087 c 


470 d 




15.40 
15.42 
15.44 

15.36 
15.38 
15.40 

15.42 
15.44 
15.46 

15.36 
15.38 
15.40 

15.42 
15.44 
15.46 


3.05 
2.51 
2.99 

2.68 
2.38 
2.32 

2.36 
2.52 
2.70 

2.90 
2.52 
2.41 

2.32 
2.41 
2.43 


4.37 
4.45 
4.54 

4.47 
4.40 
4.40 

4.47 
4.56 
4.67 

4.54 
4.46 
4.46 

4.54 
4.62 

4.72 


2.97 
3.00 
3.04 

3.01 
2.99 
2.99 

3.01 
3.05 
3.09 

3.04 
3.01 
3.01 

3.04 
3.07 
3.11 


0.106 
.092 
.089 

.104 
.092 
.086 

.079 
.077 
.078 

.094 
.089 
.092 

.095 
.098 
.106 



T 

deg°K 


d c , 
g mol/1 


(A p , %) • 100 


AC,, 

J/mol deg 


88 a 


781 b 


2087 c 


470 d 


32.982 
32.982 
32.982 


15.40 
15.44 
15.48 


2.51 
2.60 
3.11 


4.41 

4.57 
4.78 


2.99 
3.05 
3.13 


0.086 
.076 
.083 


32.984 
32.984 
P2.984 
32.984 


15.43 
15.44 
15.45 
15.46 


2.59 
2.67 

2.74 
2.88 


4.51 
4.56 
4.61 
4.67 


3.03 
3.05 
3.07 
3.09 


.077 
.075 
.076 

.077 


32.986 
32.986 
32.986 


15.40 
15.44 
15.48 


2.55 
2.69 
3.24 


4.42 
4.56 
4.78 


2.99 
3.05 
3.13 


.086 
.076 
.082 


32.990 
32.990 
32.990 


15.40 . 

15.44 

15.48 


2.74 ' 

2.88 

3.50 


4.43 
4.56 
4.78 


3.00 
3.05 
3.13 


.088 
.079 
.079 



a 88 Points, 32 ^ r i 
b 781 Points, 17 *S 
c 2087 Points, 17* 
d 470 Values of C v 



1 *S 35 °K, 8.5 *S d *s 22.5 g mol/liter. 
7^40, 1.0^^^36.0. 
= r^ 100, 1.0 *£</*£ 43.0. 
, [3] and table 5. 



Table 5. Calculated specific heats from eq (4) with n = 8, d a 15.44 g mol/l*!" 



SPECIFIC HEAT AT CONSTANT V0LUMF. JOULE/GMOL DFG. 
TEMPERATURE. DE6.K 



DENSTY 


30. 


Q0 


31. 


00 


32. 


00 


33. 


00 


34. 


00 


35. 


00 


36. 


00 


GMOL/L 


CALCD 


REFER 


CALCO 


REFER 


CALCD 


REFER 


CALCD 


REFER 


CALCn 


REFER 


CALCP 


REFER 


CALCP 


REFER 


0.000 


12.47 


12.47 


12.47 


12.47 


12.47 


12.47 


12.47 


12.47 


12.47 


12.47 


12.47 


12.47 


12.48 


12.48 


1.000 


0.35 


0.31 


0.33 


0.29 


O.30 


0.28 


0.28 


0.27 


0.26 


0.25 


0,24 


0.24 


0.23 


0.22 


2.000 


0.75 


0.70 


0.69 


0.66 


0.6* 


0.62 


0.59 


0.58 


0.55 


0,53 


0.50 


0.49 


0.47 


0.45 


3.000 


1.18 


1.17 


1.10 


1.09 


1.02 


1.01 


0.93 


0.95 


0.86 


0.85 


0.78 


0.76 


0.72 


0.68 


4.000 


1.65 


1.73 


1.54 


1.59 


1.43 


1.47 


1.31 


1.37 


1.19 


1.20 


1.0* 


1.06 


0.9« 


0.93 


5.000 


2.15 


2.37 


2.0? 


2.16 


1.89 


1.98 


1.72 


1.85 


1.55 


1.58 


1.3P 


1.37 


1.24 


1.19 


6.000 


0.00 


0.00 


2.56 


2.80 


2.41 


2.55 


2.19 


2.37 


1.94 


2.00 


1.70 


1.70 


1.50 


1.46 


7.000 


0.00 


0.00 


0.00 


COO 


3.01 


3.19 


2.74 


2.95 


2.39 


2.45 


2.04 


2.06 


1.76 


1.74 


8.000 


0.00 


0.00 


0.00 


0.00 


3.72 


3.70 


3.38 


3.60 


2.89 


?.96 


2.40 


2.35 


2.01 


2.02 


9.000 


0.00 


0.00 


0.00 


0.00 


O.OO 


0.00 


4.11 


4.33 


3.45 


3.44 


2.7« 


2.66 


2.26 


2.10 


10.000 


0.00 


O.PO 


0.00 


0.00 


o.oo 


0.00 


4.91 


4.89 


4.05 


3.80 


3.16 


2.99 


?.48 


2.55 


11.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


5.74 


5.50 


4.65 


4.37 


3.53 


3.33 


2.68 


2.76 


12.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


6.48 


6.16 


5.19 


4. 85 


3*84 


3.63 


2.84 


2.92 


13.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


7.03 


6.81 


5.58 


5.?8 


4.06 


3.82 


2.93 


3.01 


14.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


7.27 


7.26 


5.74 


5.54 


4.14 


3.88 


?.95 


3.02 


15.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


7.15 


7.36 


5.64 


5.55 


4,05 


3.80 


2.87 


2.95 


16.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


6.66 


7.16 


5.27 


5.28 


3.81 


3.32 


2.72 


?.62 


17.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


5.86 


6.36 


4.68 


4.71 


3.43 


2.99 


2*50 


2.41 


18.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


4.88 


5.23 


3.95 


3.94 


2.9* 


2.59 


2.24 


2.17 


19.000 


0.00 


0.00 


0.00 


0.00 


o.oo 


0.00 


3.86 


3.98 


3.20 


3.09 


?.50 


2.18 


1.97 


1.89 


20.000 


o.oo 


0.00 


0.00 


0.00 


o.oo 


0.00 


2.91 


2.79 


2.50 


2.29 


2.05 


1.79 


1.71 


1.61 


21.000 


0.00 


0.00 


o.oo 


o.oo 


0*00 


0*00 


2.13 


1.83 


1.91 


1 .63 


1.67 


1.46 


1.47 


1.35 


22.000 


0.00 


0.00 


0.00 


o.oo 


o.oo 


o.oo 


1.54 


1.42 


1.46 


1.31 


1.37 


1.23 


1.28 


1.15 


23.000 


0.00 


0.00 


o.oc 


0*00 


-0*07 


-0.05 


1.12 


1.17 


1.14 


1.10 


1«14 


1.05 


1»1? 


1.00 


24.000 


0.00 


0.00 


o.oo 


0*00 


-0.38 


-0.29 


0.85 


0.97 


0.93 


0.93 


0»98 


0.90 


1*00 


0.87 


25.000 


0.00 


0.00 


0.00 


0.00 


-0.58 


-0.48 


0.68 


0.81 


0.79 


0.79 


0.87 


0.79 


0.92 


0.78 


26.000 


0.00 


0.00 


-0.13 


-0.16 


-0.70 


-0.63 


0.58 


0.69 


0..70 


0.69 


0.80 


0.70 


0.86 


0.71 


27.000 


-0.03 


-0.04 


-0.22 


-0.29 


-0.77 


-0.74 


0.52 


0.60 


0.65 


0.6? 


0.75 


0.65 


0.8? 


0.67 


28.000 


-0.11 


-0.15 


-0.28 


-0.39 


-0.62 


-0.81 


0.49 


0.54 


0.6? 


0.5* 


0.73 


0.6? 


0.81 


0.66 


29.000 


-0.17 


-0.23 


-0.32 


-0.44 


-0.84 


-0.85 


0,47 


0.52 


0.62 


0.57 


0.73 


0.62 


0.81 


0.66 


30.000 


-0.20 


-0.27 


-0.34 


-0.47 


-0.84 


-0.86 


0.48 


0.52 


0.63 


0.58 


0.75 


0.64 


0.83 


0.69 


31.000 


-0.21 


-0.28 


-0.33 


-0.47 


-0.83 


-0.85 


0.49 


0.53 


0.65 


0.60 


0.77 


0.67 


0.86 


0.72 


32.000 


-0.20 


-0.27 


-0.32 


-0.45 


-0.81 


-0.83 


0.52 


0.57 


0.68 


0.64 


0.81 


0.71 


0.90 


0.77 


33.000 


-0.19 


-0.25 


-0.30 


-0.42 


-0.78 


-0.80 


0.56 


0.60 


0.72 


0.68 


0.86 


0.76 


0.95 


0.82 


34.000 


-0.16 


-0.23 


-0.27 


-0.39 


-0.75 


-0.76 


0.60 


0.65 


0.77 


0.73 


0.91 


0.81 


1.00 


0.87 


35.000 


-0.14 


-0.20 


-0.24 


-0.36 


-0.71 


-0.72 


0.64 


0.69 


0.81 


0.77 


6.95 


0.85 


1.06 


0.92 


36.000 


-0.12 


-0.18 


-0.21 


-0.33 


-0.68 


-0.69 


0.68 


0.72 


0.86 


0.81 


1.00 


0.90 


l.il 


0.97 


37.000 


-0.10 


-0.17 


-0.19 


-0.31 


-0.65 


-0.67 


0.71 


0.75 


0.90 


0.85 


1.05 


0.94 


1.16 


1.02 


38.000 


-0.09 


-0.16 


-0.17 


-0.30 


-0.63 


-0.64 


0.74 


0.78 


0.93 


0.88 


1.09 


0.98 


l.?0 


1.06 


39.000 


-0.09 


-0.16 


-0.16 


-0.29 


-0.61 


-0.63 


0.76 


0.80 


0.96 


0.91 


1.1? 


1.01 


l.?4 


1. 10 


40.000 


-0.09 


-0.16 


-0.16 


-0.29 


-0.60 


-0.62 


0.78 


0.82 


0.98 


0,93 


1*15 


1.04 


l.?7 


1.14 


41.000 


-0.11 


-0.16 


-0.16 


-C29 


-0.60 


-0.62 


0.79 


0.83 


0.O9 


0.94 


1.17 


1.06 


1.30 


1.15 


42.000 


-0.12 


-0.18 


-0.18 


-0.30 


-0.61 


-0.62 


0.78 


0.82 


1.00 


0.94 


1.18 


0.00 


1.31 


0.00 
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Table 5. (continued) 



DENSTY 

GMOL/L 

0.000 

1.000 
2.000 
3.000 
4.000 
5.000 
6.000 
7.000 
8.000 
9.000 
10.000 
11.000 
12.000 
13.000 
14.000 
15.000 
16.000 
17.000 
18.000 
19.000 
20.000 
21.000 
22.000 
23.000 
24.000 
25.000 
26.000 
27.000 
2B.000 
29.000 
30.000 
31.000 
32.000 
33.000 
34.000 
35.000 
36.000 
37.000 
38.000 
39.000 
40.000 
41.000 
42.000 



40.00 
CALCD REFFR 
12.49 12.49 



SPECIFIC HEAT AT CONSTANT VOLUMF. JOULF/GMOI 
TEMPERATURE. DFG.K 
50.00 60.00 70.00 RO.OO 
CALCD REFER CALCD REFER CALCD REFER CALCD RPFFR 
12.64 12.64 13.09 n.OR 13.99 13.99 15.31 15.31 



0.18 
0.36 
0.55 
0.74 
0.92 
1.09 
l.?4 
1.38 
1.49 
1.58 
l.ft3 
1.67 
1.67 
1.66 
1.62 
1.56 
1.49 
l.*l 
1.32 
1.24 
1.16 
1.08 
1.02 
0.97 
0.94 
0.92 
0.91 
0.93 
0.95 
0.99 
1.04 
1.10 
1.16 
1.23 
1.30 
1.36 
1.43 
1.49 
1.54 
1.59 
1.64 
1.67 



0.18 
0.35 
0.51 
0.66 
0.R0 
0.93 
1.05 
1.21 
1.38 
1.51 
1.62 
1.69 
1.73 
1.72 
1.68 
1.62 
1.52 
1.40 
1.28 
1.16 
1.05 
0.R6 
0.88 
0.82 
0.79 
0.77 
0.77 
0.78 
0.82 
0.86 
0.91 
0.97 
1.04 
1.10 
1.16 
1.23 
1.29 
1.35 
1.42 
0.00 
0.00 
0.00 



0.12 
0.24 
0.35 
0*46 
0.56 
0.66 
0.75 
0.83 
0.89 
0.9S 
0.99 
1.0? 
1.04 
1.05 
1.05 
1.04 
1.03 
1.02 
1.01 
1.01 
1.00 
LOO 
1.01 
1.03 
1.05 
1.0" 
l.U 
1.19 
1.25 
1.3? 
1.39 
1.47 
1.56 
1.64 
1.73 
1.82 
1.90 
1.99 
2,07 
2.15 
2.22 
2.29 



0.12 
0.24 
0.35 
0*44 
0.53 
0.61 
0.68 
0.75 
0.81 
0.86 
0.90 
0.92 
0.94 
0.95 
0.9* 
0.95 
0.95 
0.94 
0.93 
0.93 
0.93 
0.93 
0.94 
0.96 
0.99 
1.03 
1.07 
1.13 
1.19 
1.26 
1.33 
1.41 
1.49 
1.58 
1.66 
1.75 
1.85 
1.9ft 
0.00 
0.00 
0.00 
0.00 



0.09 
0.18 
0.?6 
0.34 
0.42 
0.49 
0.56 
0.62 
0.67 
0.72 
0.76 
0.79 
0.82 
0.85 
0.87 
0.90 
0.92 
0.94 
0.96 
0.99 
1.02 
1.05 
1.09 
1.14 
1.19 
1.25 
1.31 
1.38 
1.46 
1.54 
1.63 
1.72 
1.81 
1.91 
2.00 
2.10 
2.20 
2.29 
?.39 
2.48 
?.57 
2.65 



0.09 
0.18 
0.27 
0.35 
0.43 
0.51 
0.58 
0.63 
0.67 
0.70 
0.73 
0.75 
0.78 
0.80 
0.82 
0.84 
0.87 
0.90 
0.93 
0.96 
1.00 
1.0* 
1.08 
1.13 
1.18 
1.24 
1.30 
1.37 
1.45 
1.53 
1.61 
1.70 
1.79 
1.89 
2.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



0.0ft 
0.15 
0.22 
0.29 
0.35 
0.41 
0.47 
0.52 
0.57 
0.62 
0.66 
0.70 
0.74 
0.78 
0.81 
0.85 
0.89 
0.93 
0.97 
1.01 
1.06 
1.11 
1.16 
1 .22 
1.29 
1.36 
1.43 
1.51 
1.60 
1.6fl 
1.78 
1.87 
1.97 
2.07 
2.17 
2.27 
2.37 
2.47 
2.57 
2.67 
2.77 
2.86 



0.06 
0.13 
0.20 
0.27 
0.33 
0.39 
0.45 
0.50 
0.54 
0.57 
0.61 
0.64 
0.67 
0.71 
0.74 
0.78 
0.82 
0.86 
0.90 
0.94 
0.99 
1.05 
1.10 
1.16 
1.23 
1 .30 
1.37 
1.45 
1.53 
1.62 
1.71 
1.81 
1.92 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



0.07 
0.13 
0.20 
0.?6 
0.3? 
0.37 
0.4? 
0.47 
0.5? 
0.57 
0.6? 
0.66 
0.71 
0.75 
0.79 
0.84 
O.ftft 
0.93 
0.9P 
1.04 
1.09 
1.15 
1.?? 
1.28 
1.35 
1.43 
1.51 
1.59 
1.68 
1.77 
1.86 
1.96 
2.06 
2.16 
2.?6 
2.36 
2.47 
2.57 
2.67 
?.7fl 
?.8fl 



0.06 

o.n 

0.17 
0.22 
0.27 
0,3? 
0.37 
0.4? 
0.47 
0.5? 
0.57 
0.61 
0.66 
0.70 
0.74 
0.79 
0.P3 
0.87 
0.9? 
0.97 
1.0? 
1.0" 
1.15 
1.21 
1.28 
1.36 
1.44 
1.5? 
1.61 
1.70 
1.81 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



2.9fl O.OO 



90.00 
CALCD REFER 
16.94 16.94 



0.05 
0.10 
0.14 
0.18 
0.23 
0.27 
0.31 
0.37 
0.43 
0.49 
0.56 
0.61 
0.67 
0.72 
0.77 
0.82 
0.86 
0.91 
0.96 
1.01 
1.06 
1.13 
1.19 
1.26 
1.34 
1.4? 
1.51 
1.60 
1.69 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



loo.oo 

CALCD REFFR 
18.69 18. 69 



0.06 
0.1? 
0.17 
0.23 
0.28 
0.33 
0.38 
0.43 
0.48 
0.53 
0.5R 
0.63 
0.68 
0.73 
0.78 
0.83 
0.89 
0.95 
1.01 
1.07 
1.13 
l.?0 
1.27 
1^5 
1.4? 
1.50 
1.58 
1.67 
1.76 
1.85 
1 .94 
2.04 
2.14 
2.?4 
2.34 
2.44 
2.55 
2.65 
2.75 
2.86 
2.96 
3.06 



0.04 
0.08 
0.12 
0.15 
0.19 
0.22 
0.?5 
0.32 
0.40 
0.48 
0.56 
0.63 
0.69 
0.76 
0.81 
0,«6 
0.91 
0.95 
1,00 
1.05 
1.11 
1.17 
l.?4 
1.32 
1.40 
1.49 
1.58 
0.00 
0.00 
0.00 
0.00 
0.00 
O.OO 
0.00 
O.OO 
O.OO 
0.00 
0.00 
0.00 
O.OO 
0.00 
0.00 



FOR 470 POINTS, MEAN 
t 



0.077 J/GMOL DFG. 

At T<T the reference states for liquid are the saturated liquid (not ideal gas) with 
properties given by table 6. 



TABLE 6. Properties of saturated liquid used in the specific heat 
calculations 



T, °K 


d, 


g mol/liter [15] 


C w ,J/gmo] de-[3| 


30.0 




26.753 


13.07 


31.0 




25.103 


13.33 


32.0 




22.825 


13.79 



near d= 14.0 g mol/liter may be an imperfection, even 
though it is consistent with the estimated sign inver- 
sion of (d 2 PldT 2 ) near 13.0 g mol/liter at 33 °K (fig. 3). 
This could be the subject of another, lengthy investi- 
gation, which might include adjustment of data for 
gravitational effects, nonanalytic functions, etc. [8J. 
Further to validate the behavior of eq (4), we have 
made the following examinations. (1) A table of rela- 
tive deviations, similar to table 3 in this report, shows 
that they are comparable in the critical region with 
those elsewhere, and that many near the saturation 
envelope at T < T c have been appreciably diminished 
relative to table 3. (2) A table of pressures and deriva- 
tives dP I dp, d 2 P/dp 2 on the critical isotherm shows 
that dP/dp^0, and that d 2 P/dp 2 § as p§ p c . Dif- 



ferences in () 2 Pj()p 2 give the estimate d :i P/dp :i = 9.0 
atm at p c (where p = d/dt). (3) A table of calculated 
pressures at all p, T arguments of the interpolated 
data [5 1 shows monotonic increase of P with either 
p or T in single-phase regions. (4) A table of calcu- 
lated densities of the saturated vapor and liquid phases, 
and of the freezing liquid phase, shows acceptable 
agreement with previously derived results [15j. 

In conclusion, we mention a further examination of 
eq (4-a). Effects of varying T c (P c obtained from [13]) 
are in table 4(B). We do not feel that these results 
justify a change of T c from [12J. Constants corre- 
sponding to the minimum deviation, AC r = 0.075 J/g 
mol deg, nevertheless are given in table 7 because the 
critical point constants are close to those determined 
initially by Hoge and Lassiter, namely T c = 32.984 
±0.02 °K, d c = 15.27 ±0.30 g mol/liter, P c = 12.77 
±0.04 atm [14]. 



The National Aeronautics and Space Administra- 
tion has supported this work from the beginning. 
Mr. William J. Hall generously helped the author with 
initial computer programs, and contributed the essen- 
tial iteration routine and the high-speed, least-squares 
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Table 7. Constants for alternate T c 
Equation (4— a) 



Constant 


Value 


n 


8 


T c 




°K 


32.984 


dc 

Pc [13]... 

c 8 ... 




g mol/liter 

atm 


15.44 

12.774 257 
-0.7702 5576 


D x 


.0004 4993 


/),..., 


.01614536 


Di 


- .0102 2498 


D 4 


7.22 


£,... 


3.32 







routine in machine language. Mrs. Levelt Sengers 
kindly suggested improvements in the manuscript 
and encouraged us to seek a form which would be non- 
analytic at the critical point. 

6. Appendixes 

6.1. Method Used To Find Data for A(p), B(p), C(p) 

With reference to eq (3) and determination of coef- 
ficients on the right side by least squares on isochores, 
we used a method whereby deviations introduced by 
the arbitrary form of the first coefficient, A(p), will ap- 
pear in data for the second coefficient, B(p), and so 
forth. After adjusting the constants in A(p) we used 
its representation on the left side, and redetermined 
data for B and for C by least squares on isochores in 
the transposed form, 

{log Y— A ca{c • p - x X ' 4 }=B - p - x~ l + C • x~ 3 9 

and so forth. 

6.2. The Requirement for High Precision in A(p), Etc. 
The explicit form of eq (3) for pressure, 

P=-U(p)+F(p,T) 

shows that P is obtained as a difference of two prin- 
cipal terms. At high densities the value of II may 
exceed 700 atm. At high densities and low tempera- 
tures the value of P may be as low as 0.0695 atm (triple 



point). Values for the parameters A(p), B(p), etc., 
therefore may be required with a precision ten thou- 
sand-fold greater than desired in the calculated 
pressure. 
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